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Abstract 

A three-parameter family B = B(a,b,c ) of weighted Hankel matrices is intro¬ 
duced with the entries 

_ T(j + k + a) /r(j + b)T(j + c)T(k + b)T(k + c) 

3 ' k T(j + k + b + c) y T(j+ a)j\T(k + a)k\ 

j, k £ Z + , supposing a , b, c are positive and a < b + c, b < a + c, c < a + b. 
The famous Hilbert matrix is included as a particular case. The direct sum 
B(a, b, c) © B{a + 1, b + 1, c) is shown to commute with a discrete analog of the 
dilatation operator. It follows that there exists a three-parameter family of real 
symmetric Jacobi matrices, T(a,b, c), commuting with B(a,b,c). The orthog¬ 
onal polynomials associated with T(o, b, c) turn out to be the continuous dual 
Hahn polynomials. Consequently, a unitary mapping U diagonalizing T(a,b,c ) 
can be constructed explicitly. At the same time, U diagonalizes B(a,b,c) and 
the spectrum of this matrix operator is shown to be purely absolutely contin¬ 
uous and filling the interval [0 , M(a,b,c)\ where M(a,b,c ) is known explicitly. 
If the assumption c < a + b is relaxed while the remaining inequalities on a, b, 
c are all supposed to be valid, the spectrum contains also a finite discrete part 
lying above the threshold M(a, b, c ). Again, all eigenvalues and eigenvectors are 
described explicitly. 
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1 Introduction 


An integral operator K on L 2 ((0, oo), dr) whose integral kernel K,(x,y) is real, sym¬ 
metric and homogeneous of degree —1 on the first quadrant and such that 

|/C(t, 1)| f -1//2 df < oo 

is bounded and explicitly diagonalizable by the Mellin integral transform. In more 
detail, K is unitarily equivalent to the multiplication operator by the function 

f-OO p 

g(0= / /C(t, 1) r 1/2 "*dt = / e -^ x JC{e x / 2 ,e- x/2 )dx 

Jo Jr 

acting on L 2 (M, d£). This feature can readily be understood if the symmetry properties 
of such an integral operator are examined. K commutes with the one-parameter 
unitary group of dilatation transformations on the positive half-line generated by the 
skew-symmetric differential operator D = xd/dx + 1/2. The kernel of the Mellin 
integral transform is in fact nothing but a family of generalized eigenfunctions of D. 
All what has been said above is applicable to the integral kernel 



£e(x,y) 


(. xy ) e / 2 
(x + y) l+1 


( 1 ) 


depending on a real parameter £ provided £ > — 1. The corresponding integral operator 
Kp is unitarily equivalent to the multiplication operator by the function 


9(0= / t ( *“ 1)/2 +*(t + 1) 

Jo 


-i-i 


dt = 


T(£ + 1) 


ri -(* + !) + >£ 


( 2 ) 


In no way it is straightforward to find an authentic discrete analog of the integral 
kernel (IT]) . Of course, given a homogenous kernel of degree —1 one can always restrict 
the kernel to the discrete set (0 + Z + ) x (9 + Z + ), for some 9 > 0 and with standing 
for nonnegative integers, obtaining this way a semi-infinite matrix. For example, using 
the kernel (HD, with t = 0, we get the (generalized) Hilbert matrix. Matrices of this 
type have been explored, for instance, in [7j. But as emphasized in the introduction of 
the cited paper, it appears that there may be inconveniences in applying to matrices 
some methods originally invented for integral operators. In particular, let us note that 
no obvious discrete analog of the Mellin integral transform yielding a diagonalization of 
such matrix operators is at our disposal. Nevertheless, comparison of a matrix operator 
obtained by discretizing an integral kernel with the integral operator in question can 
provide, under certain additional assumptions, quite a useful information including 
the precise value of the norm, see Clip. 1 in H. 

The present paper actually aims to describe matrix operators whose character and 
properties remarkably closely resemble those of the integral operators K(. Notably, 
despite of nontrivial form of such matrices all spectral properties are known in full 
detail. Our approach is based on a construction of a three-parameter family of matrices 
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B = B(a, b, c ) with the desired structure of weighted Hankcl matrices and enjoying an 
appropriate symmetry. By the symmetry we mean a discrete analog D of the generator 
of dilatation transformations which is commuting with B. D can be regarded as a 
first order difference operator. Then D 2 commutes with B as well and, since this is 
a second order difference operator, as such can be related to a Jacobi (tridiagonal) 
matrix. Owing to the discrete symmetry D, an explicit diagonalization of B is possible. 
To this end, one has to rely on the theory of orthogonal polynomials in place of the 
Mellin transform. This is only roughly the main idea which is made precise and fully 
developed in the remainder of the paper. 

The famous generalized Hilbert matrix is included in the family B(a,b,c ) as a 
one-parameter subfamily. The constructed diagonalization procedure if applied to 
Hilbert’s matrix reproduces, in an alternative way, a complete solution to the spectral 
problem due to Rosenblum ra- Considering only the original Hilbert matrix, i.e. not 
treating the whole one-parameter family, it has been observed earlier by Otte that a 
commuting Jacobi matrix can be used for a diagonalization. This result exists in the 
form of a conference presentation which is currently available from author’s web page 
[ID] . But the first one who has discovered that the generalized Hilbert matrix has a 
tridiagonal matrix in its commutant seems to be Griinbaum in [6]. As noted in the 
cited paper, the class of Hankel matrices with this property is very limited. Saying 
this, the author was guided by his forgoing studies of a similar problem for the class 
of Toeplitz matrices |5J. 


2 A family of weighted Hankel matrices with a dis¬ 
crete symmetry 

Consider the three-parameter family of semi-infinite symmetric matrices B = B (a, b, c), 


Bj,k 


T(j + k + a ) [rJj + b)T(j + c)T(k + b)T(k + c ) 

T{j + k + b + c) y T(j + a) j\ T(k + a) k\ 


(3) 


The parameters a, b and c should be restricted to a range for which the matrices 
B(a,b,c ) are real and hence Hermitian. Throughout the paper we shall assume that 
the parameters a, b , c are all positive. The entries of B have the structure Bj ik = 
w(j)w(k)h(j + k ) and therefore the matrices may be classified as weighted Hankel 
matrices, see Clip. 6 §8 in DU. 

Note that, by the Stirling formula, the leading asymptotic term of the matrix entry 
Bj k for both j and k large is 


Bj t k 




(j'/ c )(6+c-a- 1 )/ 2 

(j + k) b+c ~ a 


This may be compared with the integral kernel ([!]), with £ = b + c — a — 1, suggesting 
that we should suppose 

b + c — a > 0. (4) 
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The matrix entries B jk take particularly simple form if we put a = 9, b = 6 + £, 
c — 1, with £ G Z + . Then 


p _ \/ (j + 9 )e(k + 0 )^ 

(j + fc + 0)m 

where (a)i is the usual Pochhammer symbol. In particular, for l = 0 we obtain the 
generalized Hilbert matrix, Bj tk = 1 /(j + k + #); for a short account on its history see 
Clip. IX in [Sj. 

Moreover, letting £ = 1 we get a modification of the so called Bergman-Hilbert 
matrix, 

y /(J + B)(k + e) 

(j + k + e)(j + k + e + iy 

studied in m- 

Put A = B(a, b, c) © B(a + 1, b + 1, c). A can be identified with the semi-infinite 
matrix with the entries 

A'2j.2k — B(a , b, c)j tk , A2j+i,2k+i — B(a + 1,6+1, c)j jk , j,k = 0,1, 2,..., 

and A, j k — 0 if j, k are of different parity. This identification corresponds to the direct 
sum 

£ 2 (Z + ) = £ 2 (2Z + )@£ 2 (2Z + + 1). (5) 

(£ 2 (2Z + ) is spanned by k] k G Z + j, £ 2 (2Z + + 1) is spanned by {e 2 fc+i; k G Z + j 
where {e k ', k G Z + j is the standard basis in £ 2 (Z + )). Furthermore, let us introduce 
another three-parameter family of matrices, D = D(a,b,c), such that 

Dj,j+i = ~Dj+i,j = d(j), D jtk = 0 otherwise, 


where 


d(2j) = y/(j + a)(j + 6), d(2j + 1) = sj (j + l)(j + c), j,k = 0,1, 2,.... 


D can be regarded as a discrete analog of the dilatation operator. 

With respect to the decomposition (151) . the matrices A and D can be written in 
the blockwise form 


( B(a,b,c) 0 \ n _ ( 0 C \ 

V 0 B{a + l,b+l,c) )' \ -C T 0 /’ 


( 6 ) 


where C is another semi-infinite matrix with the entries 


Cjj = d(2j), Cj + \j = —d{2j + 1), Cj >k = 0 otherwise. 
Lemma 1. A and D commute. 
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Proof. The proof can be carried out relying just on a straightforward evaluation of 
the matrix entries of the commutator. One has (with d(— 1) := 0) 

(AD - DA) jtk = d(k - l)A jt k-i - d(k)A j>k+1 + d(j - 1 )Aj_ hk - d(j)A j+hk . 

Obviously, (AD — DA)j ik vanishes whenever j and k are of the same parity. Without 
loss of generality it suffices to consider the situation when j is even and k is odd. Let 
us suppose j, k G Z + and evaluate the expression 

(AD - DA) 2 j, 2 k+i = d(2k)B(a, 6, c) hk - d(2k + 1 )B(a, b, c ) jik+ i 

+ d(2j — 1 )B(a + 1, b + 1, c)j—\ jk — d(2j)B(a + 1,6+1, c)j jk . 

A direct computation actually shows that the expression equals zero. 

One may prefer, however, to proceed another way revealing the algebraic structure 
behind the identity. In view of 0, AD — DA = 0 is equivalent to 

B(a, 6, c) C = CB(a + 1, b + 1, c), (7) 

i.e. C intertwines the operators B(a,b,c ) and B(a + 1,6 + l,c). In order to take 
into account the structure of B — B(a,b,c ) let us write B = WHW where W is the 
diagonal matrix with the diagonal entries 


W jd = 


T (j + b)T(j + c) 

r (j + a)j\ 


and H is the Hankel matrix, H jk — h(j + k), where 


h(z) 


T(z + a) 
T(z + 6 + c) ’ 


Similarly, B(a + 1,6 + l,c) = WHW , with analogous expressions obtained just by 
shifting the parameters a and 6 by 1. Note that h(z) = h(z + 1). Then (0 can be 
rewritten as 

HV = VH , with V = WCW~\ V = W~ X CW. 

One readily finds that 


Vjj =j + a, V j+ ij = - j - c, V jd = j + 6, V j+ ij = -j - 1, Vj- fc = V j)k = 0 otherwise. 
Consequently, 0 means that 

h(j + k)(j + k + a) = h(j + k) (j + k + 6 + c) 


for all indices j, k. But this is actually so since 

h(z) = h(z + 1) = z ^ h ° +c h(z), 

thus concluding the proof. □ 
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Corollary 2. Let T = T(a, b, c ) be the three-parameter family consisting of symmetric 
Jacobi matrices with the entries 


T j,j = 3 U + c - 1 ) + (j + a)(j + b), T jd+ 1 = T j+1J = -y/(j + + a )U + b )U + c )> 

Tj : k = 0 otherwise. Then the matrices B(a,b,c ) and T(a,b,c ) commute. 

Proof. Clearly, A and D 2 commute and therefore, in view of (|6j) . B and T = CC T 
commute. T is a Jacobi matrix with the entries 

T jtj = d(2j - l) 2 + d{2jf — j (j + c — 1) + (j + a)(j + b), 

Tjj+i = Tj+i,j — —d(2j)d(2j + l) = —a / (j + 1 ){j + a)(j + b)(j + c), 

and Tj t k = 0 otherwise. □ 


3 Diagonalization and the spectral properties of 

B(a, b , c) 


3.1 The associated orthogonal polynomials 

The monic orthogonal polynomials associated with a Jacobi matrix T are defined by 
the recurrence: P-i(x) = 0, Pq(x) = 1, and 

p j+i(x) = {x- T j+1J+1 ) Pj{x) - (T jd+1 f Pj-i(x). 

In our case, it is convenient to modify the matrix T by adding a multiple of the unit 
operator. The redefined diagonal entries read 

Tjj = dU + c- 1 ) + (J + a)(j + b) - ^(a + b-c) 2 

= ~j (j + 1 ) + (j ~ b + d)(j - c + d) + (j - a + d)(j - b + d ) 

+ (j ~ a + d) ( j - c + d) 


where d = (a + b+c)/2. The associated monic orthogonal polynomials, P n {x), coincide 
with the continuous dual Hahn polynomials; see Eq. 9.3.5 in [9]. We have (see also 
0 Eq. 9.3.1]) 


2 b + c — a a + c — b a + b — c 


P n (x 2 ) = (-1 ) n S n \x*- 


. b+c—a b+c—a 
= (-1) (6) n (c) n 3 F 2 -n,--- htx, --- ix]b,c] 1 


The measure of orthogonality for the continuous dual Hahn polynomials is known 
explicitly 0 Eqs. 9.3.2, 9.3.3]. Put 


p(x) = 


xsinh( 27 ra;) 

7T 2 r(a)r(5)r(c) 


b + c — a \ ^(a + c 

---h ix ) T f 


b . \^f a + b — c 
-h ix 1 ---h ix 
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Then, for a, b, c positive such that b + c > a, a + c > b and a T b > c, we have the 
orthogonality relation 

).P n (x )p(x) (\x (o7) n (JX) n (^C) n Til b m n . 

More generally, suppose for definiteness that 0 < b < c and, in agreement with (dJ) , 
0 < a < b + c. Then obviously a + c — b > 0 but it may happen that a T b — c < 0 
(then necessarily b < c). In that case the orthogonality relation should be modified 
by adding a discrete part. So suppose 

0 < b < c and 0 < a < c — b (8) 

(then a < c — b < b + c). Put (|br] meaning the ceiling of x G M) 

N(cl , 6, c) = |~(c — ci — 6)/2~| — 1. (9) 



One has 


P m (x 2 )P n (x 2 )p(x ) dx + 


T(c — a)T(c — b) 
Y{c)Y{c-a-b) 


N(a,b,c) 

x £ 


k =0 



(—l) fc (a + b - c) k {a)k{b)k 

(a — c + 1 ) k (b - c+ 1 ) k k\ 


X Pm 


a + b — 
2 


c 




(®)n(^)n(c)n Til Sm,n- 


( 10 ) 


3.2 Diagonalization of B{a, 6, c) 

Assuming ((3]) let 


Pjx 2 ) = S„ 


x 


b + c — a a + c — b a + b — c 


\J (a)n(b) n (c) n n\. ( 11 ) 


The polynomials P n (x 2 ) are normalized and ( K P 0 (x 2 ), Pi(x 2 ), P 2 (a; 2 ), • • •) is a formal 
eigenvector of the matrix operator T(a,b,c ) corresponding to the eigenvalue x 2 . Let 
us introduce the unitary transform 


U : £ 2 (Z + ) — > L 2 (M(a, b, c), dp) : e n —> P n (x 2 ), n G Z + , (12) 

where 

M(a , b , c) = (0, Too) U {A fc ; k = 0,1,..., N(a, b , c)}, X k = i H (13) 


(the discrete part occurs if and only if a T b — c < 0), dp(x) = p(x)dx on (0, Too) and 

_ (—l)*T(c - a)T(c - b) / 2k \ (a + b - c) k (a) k (b) k 

! k T(c)T(c — a — b ) \ a + b — cj (a — c + l) k (b — c + l) k k\ 
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Remark. Note that if 0 < k < N(a, b , c) and (JHJ) is true then /u({Afc}) > 0 as it should 
be. In fact, the signs of the numbers (a + b — c)k, (a — c + 1+ (b — c + l)fc are all equal 
to (— l) k . This is also a standard fact that (Ai(a,b,c), dp) is a probability space, as 
seen from (with m — n — 0). 

Theorem 3. The matrix operator B(a,b,c ) on £ 2 (Z + ) is unitarily equivalent to the 
multiplication operator by the function 


h{x) 


1 

r(6 + c 




(14) 


acting on L 2 (M. (a, 6, c), dp). 


Remark. Once more, this result should be compared to the integral operator K f which 
is unitarily equivalent to the multiplication operator by the function g(£) introduced 
in (J2J) and acting on L 2 (M, d£). Again, one has to put £ = b + c — a — 1. 

Proof. Directly from the construction it follows that U diagonalizes T, namely UTU -1 
equals the multiplication operator by x 2 . Since B and T commute and the spectrum 
of T is simple, UBU^ 1 is necessarily a multiplication operator, too, say by a function 
h(x)] see, for instance, Lemma 6.4 in [133 or Proposition 1.9 in Supplement 1 of |2]. 
One derives 

OO 

h(x) = h(x)P 0 (x 2 ) = UBe 0 = 57 B 0 jPj(x 2 ) 

j=o 

E r(b)r(c) _ / 9 b + c — a a + c — b a + b — c\ 

_ r(j + 6 + c) ~j\ s '0 ; -2-' —2—- -T-) ■ 


More conveniently, one can rewrite the expression in terms of the Beta function, 


K x ) = J 2 


B (b,j + c) f 2 b + c — a a + c — b a + b — c 


3=0 




SAx 2 - 


and so 

h( x ) = I (i -ty 1+b t 


1 / OO 

— 1+6 +— 1 +c 


E P f 2 b + c — a a + c — b a + b — c 


d t. 


Making use of the generating function (see Eq. 9.3.12 in [9]) 


J-71 

57 7——7 S n (x 2 ;a,f3, 7 ) = (1 - t)~ 1+lx 2 F x (a + ix,/3 + ix\a + ft;t) 


one has 


h{x)= I (l-t)~ 1+ ^ b+c ~ a) F+ ix t~ 1+c 2 F 1 i l ^P-—- + ix, —- + ix;c;t)dt. 





















Hence 


. . \ ■> 1 fb T c — o A ( cl c — b 

h{x) = ) -7 I --- 1 - ix 1 I ---h ix 

_ n ( c)»«! 


n =0 


-1 


y / l+(6+c—a)/2+ix^—1+n+c ^ 

Jo 

fb + c-a . 

= BI --- +?x,c 

fb+c—a a+c—b b+c—a 

x 2_Fl I —2-^ vx ’ —2-^—2-^ c + ^; 1 ) • 

Recalling that [TJ Eq. 15.1.20] 

2 Fi(a,/ 3 ; 7 ; 1) = - Q —^ if Re(q - a - 0) > 0 

E(t - a) r (7 - /H 

we finally obtain the desired expression. 


□ 


3.3 The spectrum of B(a,b,c ) 

The function h(x) defined in (fl4l) is bounded on the positive half-line and therefore the 
matrix B(a, b, c ) represents a bounded operator on £ 2 (Z + )- In what follows, B(a, b, c ) 
is interpreted in this manner. 


Corollary 4. The absolutely continuous part of the spectrum of B(a,b,c) is simple 
and fills the interval [0, M(a, b, c)] where 


M(a, b , c) 


1 r ( b + c ~ a \ 2 

T(b + c-a) \ 2 ) 


Proof. From Theorem [3] one infers that the absolutely continuous spectrum of B (a, b, c ) 
fills the closure of /i([0, +oo)). We have h( 0) = M(a,b,c ) and h(+o o) = 0. Moreover, 
for any u real fixed, the function f(x) = |T(ti + ia ;)| 2 is monotone decreasing on 
(0, Too). This is immediately seen from the product formula for the Gamma function 
[U Eq. 6.1.3] yielding 


1 

|T(u T ix )| 2 


(u 2 + x 2 ) e 2lu J [ 

71=1 



g —2 u/n 


where 7 is Euler’s constant. The assertion follows. 


□ 


Suppose (JHJ) and recall the notation introduced in (JHJ) , (fl3|) . Then the point spec¬ 
trum of T consists of the points X 2 , k — 0,1, ..., N(a, 6 , c), each eigenvalue is simple 
and for an eigenvector corresponding to X 2 one can choose 

OO 

v k = y ^ Pj{\) e j- 

i=o 

These vectors are also eigenvectors of B. 


9 















Corollary 5. Assuming fffif and that a + c — b, a + b — c are both nonnegative, the 
point spectrum of B(a,b,c) is empty. Assuming the point spectrum of B(a,b,c) 
equals 

spec p B(a , b, c ) = {(3 0 , ..., p N (a,b,c)} 

where 

a _ u/ \ \ _ r (6 + k)T(c - a - k) 

Pk ■— rt{A k ) - - —— - 7 -, 

1 (o + c — a) 

and it holds true that 

Po> Pl> ■■■> fi N(a,b,c) > M(a, b, c ). (15) 

In particular, all eigenvalues are simple. For an eigenvector corresponding to f3 k one 
can choose the vector v k with the components 


{p n ,Vk) P n { A fc ) 


Its norm fulfills 



3 F 2 (—n, b + k, c — a — k; b, c; 1 ). 


ikir 


r(c)r(c — a — b — k + 1) k\ 

(c — a — b — 2k)T(c — a — k)T(c — b — k)(a) k (b) k 


(16) 


Proof. Theorem [3] implies that the eigenvalues of B(a, b, c ) are exactly the values h(X k ). 
k = 0 , 1 ,..., N(a, b, c ). It remains to show (fT51) . As far as (TT5T) is concerned, one 
has 

^fb + c — a fa + b — c 
T(b + c-a)P k = T (---+)--- 



(recall that a + b — c + 2k < 0 for 0 < k < N(a , b,c)). Hence it suffices to observe that 
the function /(f) = log(T(f)) is convex on the positive half-line. It is so because 


m = m = 


oo 


E 


1 

(k + f ) 2 


where is the digamma function. It follows that, for any u > 0, the function g(t) = 
T {u — t)T(u + t ) is strictly increasing on the interval [0, w). 

Concerning (TT 6 j) . let f k G L 2 (Ai(a, b, c), dfj,(x)) be defined as follows: f k (Xj) = 5 k j 
and fk(x) = 0 on (0, oo). Then Uv k = cf k for some c G C. One immediately fords 
that 

OO 

c = (Uv„)(\ k ) = Y. hot) 2 = IMI 2 > 0- 
3=0 


On the other hand, 


v k \\ 2 = \\Uv k \\ 2 = c 2 \\f k \\ 2 = c 2 p({ X k }) = ||^||V({A fc })- 
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Whence 

1 


INI 2 


- M{^}) 

(-l) fc r(c-a)r(c-6) 

T(c)T(c-a-b) 


1 + 


2k 


(a + b - c)k(a) k {b)k 


a + b — c) (a — c + l)k(b — c + l)k k\ 


(c — a — b — 2k)T(c — a — k)T(c — b — k)(a) k (b) k 
T(c)T(c - a - b - k + 1) k\ ' 


This shows (TT 6 ]i . 


□ 


Remark 6 . Still assuming flTJ), B(a, b, c ) is a positive bounded operator on £ 2 (Z + ) and 
one has 


\\ B (a,b,c)\\ = 


T (6 + c — a) 


b + c — a 


(17) 


if a + fe — c> 0 , a + c — 6>0 (at least one of the expressions is necessarily positive) 
and 

||B(a,6,c)|| = ^ (<: ~ a) (18) 

f (b + c — a) 

if a + 6 — c<0, a + c — 6>0, and similarly if a + 6 — c > 0, a + c — fe < 0. This means 
that for every square summable real sequence {£fc}, 

OO OO / OO \ 

o )> 

j =0 k =0 \k =0 / 

and the bound is best possible. Equivalently one can also say that for any real sequence 
{£*;} and all n G Z+, 


°*EE 

j =o k =o 


T (j + k + a) 
T(j + k + b + c ) 


£j£k < \\B(a,b,c)\\ ( Y 


,fc= 0 


T(k + a) k\ 
T(k + b)T(k + c] 


& , 


with \\B(a, b, c)|| being specihed in (ITT)) . (fl 8 l) . 

The Hilbert double series inequality is a particular case for a = 6, b = 9 and c = 1 
assuming that 9 > 1/2. Explicitly, for any real square summable sequence {£*.}, 


osEE^bME* 


j =0 k =0 


j + k + 9 


K k =0 


For a = 9, b = 9, c = 1 and 0 < 9 < 1/2 one gets the inequality 


osEEt 




< 


7T 


j=0 k =0 

Again, all bounds are best possible. 


j + k + 9 sin( 7 r 6 l ) 


E« 

o 
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4 Hilbert’s matrix and the Bergman-Hilbert ma¬ 
trix 


4.1 Hilbert’s matrix 

As already remarked above, H{9) := B(9,9, 1) is the generalized Hilbert matrix, 

H{0)j,k = * + 0 , j,k = 0,1,2,.... (19) 

By our assumptions on the parameters, 6 is positive. By Corollaries [4] and [5j the abso¬ 
lutely continuous part of the spectrum is simple filling the interval [ 0 , 7 r] independently 
of 9. The point spectrum is nonempty if and only ifO<0<l/2 and if so it consists 
of the single simple eigenvalue /?o = 7r/ sin(7T0). 

Observing, however, that the dehning expression for H{9) is free of square roots, 
the range of 9 can naturally be extended to 9 G M\(—Z + ). The diagonalization 
method, as exposed in Section El can be applied to H(9) without essential modifica¬ 
tions even with this extended range. This is why we confine ourselves just to sketching 
some basic steps. 

First of all, H{9) commutes with the Jacobi matrix T{9) with the entries 


mu = 2j u+e) - 1/4+ e, T( 9 )u+i = m i+ u 


-(i + 1)0+9), 3 = 0 , 1 , 2 ,..., 


and T(9)j t k = 0 otherwise. Referring to (fTTjl . the associated normalized orthogonal 
polynomials are given by 


Pntf) = 


n\ (9), 

(0)n 


5„p 2 ;-i + e,i i 


n\ 


3 r 2 


-n, — i + 9 + ix, ~ + 9 — ix] 9,9] 1 ). 


It is useful to observe that the polynomials P n (x 2 ) can also be expressed in terms of 
the Wilson polynomials [15j . By definition, for n G Z + , 


JF n (x 2 ;a,/3,7, 6) 

(a + f})n{OL + 7 )n(cii + S) n 


4 F 3 (— n, n+a+/3+ 7 + 5 — 1 , a+ix, a—ix ; a+{3, a+ 7 , a+5 ; 1). 


By inspection of [9] Eqs. 9.1.4, 9.1.5] one finds that 



4 n (x 2 1 0 11 9 3\ 

n\ (9)'2n Wn w ; _ 4 + 2 ’ 4’ 4 + 2 1 ’ 4 ) ' 


According to (flUl) . if 9 < 1/2 then the orthogonality relation reads 



( 20 ) 


( 21 ) 
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where N(9) = [- 1/2 - ff\, 


2 


pO) 


2x tanh( 7 T.X') 

r( 0 ) 2 


r( —- + 9 + ix 


and 


A/t — i ( — — + 9 + k j, p({Afc}) 


T(l-P) 2 (l- 2 P- 2 fc) 
k\ T (2 — 29 — k) 


The sum on the LHS of (ITT!) is absent if 9 > 1 / 2 . 


Remark 7. Strictly speaking, the orthogonality relation, as described for instance in 
[9, Eqs. 9.3.2, 9.3.3], covers only the cases when 9 > 0. Nevertheless, making use of 
(15U1) and a very general complex orthogonality relation for the Wilson polynomials 
stated in ra, one can quite straightforwardly extend the desired formula to all values 
6 < 1/2, — 9 7L+. Given a, (3,^,5 G C, write for short W n (z) = W n (z; a, (3, 7 , S), 
n G Z_|_. As proved in [15 ]. 


where 


1 

2 m 



f(z)f(-z)W m (z 2 )W n {z 2 ) dz = 5 m , n MR 


f(z) 

M 


T(a - z)T{P - z)Y(j - *)r(<$ - z) 

2 r(a + P)T(a + 7)T(a + 5)T(/3 + 7 W + W7 + &) 
T(a + P + 7 + 5) 


and 

_ n! (a + /3 + 7 + 5 — 1 ) n {pt + [3) n (a + 7 )^( 0 ; + S) n ((3 + ^) n {(3 + ^) n (7 + ^)n 

(a + f3 + 7 + 5 ) 2 ^ 

The contour C is the imaginary axis deformed so as to separate the set of poles of f(z) 
from the set of poles of f(—z), assuming these two sets to be disjoint. In particular, if 
a, /3, 7 and 6 are positive, C may be taken to be the imaginary axis. The orthogonality 
measure for the Wilson polynomials is then positive and supported on the positive real 
half-line. In Section 3 of [15] this result is extended to the case when a is negative 
while a + (3, a + 7 and a + S are positive. The poles at z — ±(a + k), with k G Z + 
and a + k < 0 , then give rise to mass points of the orthogonality measure which are 
located on the negative real half-line. This procedure can readily be adapted to our 
case, with a = —1/4 + 9/2, (3 = 1/4, 7 = 1/4 + 9/2 and 5 = 3/4, finally resulting in 
the orthogonality relation (151]) . 

Relying on (151]) one can show, similarly as in Theorem [31 that H{6) is unitarily 
equivalent to the multiplication operator by the function 


cosh(7rx) 
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acting on L 2 (A4(9), d/i) where A4(9) = (0,+oo) U {A*,; k = 0,1,..., N(9)} and 
dp(x) = p(x)dx on (0, +oo). The discrete part of A4(9) occurs if and only if 9 < 1/2. 
The corresponding unitary mapping has an analogous form as that in (TT2lh From 
this explicit diagonalization one immediately deduces the full information about the 
spectral properties of H(9) thus reproducing the original result due to Rosenblum as 
stated in [12], Theorem 5]. Rosenblum’s approach was quite different than ours though 
an appropriate symmetry was heavily employed, too. Namely, it has been shown that 
H{9 ) is unitarily equivalent to an integral operator on the positive half-line such that 
there exists an explicitly diagonalizable Sturm-Liouville operator in its commutant. 

Theorem 8 . For all real 9,9 7 ^ 0, —1, —2,..., the singular continuous part of the 
spectrum of Hilbert’s matrix H (9) is empty and the absolutely continuous part is simple 
and fills the interval [ 0 , 7 r], For 9 > 1/2, the point spectrum of H(9) is empty. For 
9 < 1/2, let N(9) = [—1/2 — 9~\. Then the only possible eigenvalues of H{9) are 
it / sin( 7 r 6 ) ) and —it/ sin/yr^) whose multiplicities are respectively equal to N{9)/2 + 1 
and N{9)/2 for N(9 ) even, and they are both equal to (N(9) + l)/2 for N(9) odd. 


4.2 The Bergman-Hilbert matrix 

As another application of the general results stated in Theorem [3] and Corollaries 0] 
and [5] let us consider the so called Bergman-Hilbert matrix A with the entries 




V(j+ !)(*+!) 
(j + k + l) 2 


j,k e z+, 


which has been introduced and studied as an operator on P{fL + ) in [4,3]. It is shown 
in [3, Prop. 2] that the essential spectrum of A equals the interval [0,1]. We can make 
this analysis more complete by identifying the absolutely continuous spectrum of A. 


Proposition 9. The absolutely continuous spectrum of the Bergman-Hilbert matrix 
A, regarded as an operator on £ 2 (Z + ), is simple and fills the interval [0,1]. 

Proof. Referring to (|3]), let B = B( 1,1,2). Then 


Bj,k 


P(j + i)(k + i) 

(j + k + l)(j + k + 2)’ 


j,k e z+, 


From Corollary [4] we know that the absolutely continuous spectrum of B is simple 
filling the interval [0,1]. 

Let us show that Z := A — B is a trace class operator. We have 


Expanding 


V(j + l)(fc + l) 

3 ’ k c T + k + m 3 + k + 2y 

1 ( s+l)fk s 

(j + k + l) 2 (j + l) s+2 0 + 1) 


s+2 
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we can write 


Z = T s , where T s = (s + 1) J s i7(l) J s , J s = diag 

s=0 

where we have used the notation f[T9]i . From Theorem [8]we know that H( 1) is positive 
and so is T s . Consequently (|| • ||i standing for the trace norm), 

OO OOOO OO OO / , -1 \ -Oc 9 

s=0 s=0 j =o z j =o s=o + ' 

The space of trace class operators is complete and therefore Z is, too, a trace class 
operator. 

To conclude the proof we recall that the absolutely continuous spectrum is known 
to be invariant with respect to trace class perturbations. □ 


(j + l ) s + 3 / 2 


;jez + 
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